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Abstract The objective of the Dawn gravity investigation is to use high precision X-band 
Doppler tracking and landmark tracking from optical images to measure the gravity fields 
of Vesta and Ceres to a half-wavelength surface resolution better than 90-km and 300-km, 
respectively. Depending on the Doppler tracking assumptions, the gravity field will be de- 
termined to somewhere between harmonic degrees 15 and 25 for Vesta and about degree 
10 for Ceres. The gravity fields together with shape models determined from Dawn’s fram- 
ing camera constrain models of the interior from the core to the crust. The gravity field is 
determined jointly with the spin pole location. The second degree harmonics together with 
assumptions on obliquity or hydrostatic equilibrium may determine the moments of inertia. 

Keywords Dawn • Vesta • Ceres • Gravity • Geophysics 


1 Introduction and Background 

The Dawn mission provides the first opportunity to investigate the interior structures of two 
dwarf planet sized asteroidal bodies, namely 4 Vesta (265-km mean radius) and 1 Ceres 
(476-km mean radius). These asteroids are large, “endmember” examples of planetary em- 
bryos, with Vesta once significantly melted and volatile poor, and Ceres less processed in an 
igneous sense and with a significant volatile content. The interiors of these bodies contain 
treasure troves of information regarding the range of processes that contributed to accretion 
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of solar system objects from terrestrial planets to ice-rich outer satellites. Previous grav- 
ity investigations of other solar system objects have shown varied internal structure of the 
larger planets (e.g., Zuber et al. 1994, 2000), and the near uniform density distribution of the 
smaller asteroid 433 Eros (Konopliv et al. 2002). The primary objective of the Dawn gravity 
experiment is to measure the gravity fields of Vesta and Ceres to 90-km and 300-km half- 
wavelength resolutions, respectively. Gravity science will provide the primary information 
on the internal structures of these preserved protoplanets. 

Gravity science will be extremely challenging in the Dawn mission because of Dawn’s 
solar electric spacecraft. The frequent thrusting of the solar electric engines represents a 
non-conservative force that needs to be understood (and in the current mission plan, avoided 
for extended times) to estimate gravity. Other non-conservative forces such as solar radiation 
pressure are also increased due to the larger solar arrays required to power a solar electric 
mission. Gravity is essential to Dawn’s science objectives of characterizing the structure 
and evolution of Vesta (Zuber et al. 2011) and Ceres (McCord et al. 2011). Measurements 
of gravitational perturbations will provide constraints on the radial distribution of internal 
mass, including core size and density. Measurement of the orientation of the rotational pole 
and spin rate may constrain the polar moment of inertia from the measured obliquity (Bills 
and Nimmo 2011). Although the accuracy and amount of planned tracking data will enable 
detection of a wobble of the rotation axis, as might be excited by an impact, it is not expected 
that a detectable wobble will be present, due to the short post-collision damping time (Peale 
1976). 

Gravity, in combination with improved knowledge of the shapes of Vesta and Ceres (Ray- 
mond et al. 2011b) from observations obtained from Dawn’s framing camera (Sierks et al. 
2011), will provide insight into the distribution of crustal thickness, mechanical properties 
of the lithosphere, the possible presence of mascons and the nature of impact processes in 
general. 

Dawn’s gravity and topography will permit collectively determination of the bulk 
densities of Vesta and Ceres to better than 1%. The volume of Vesta was found from 
Hubble images with an uncertainty of about 12% (Thomas et al. 1997). The most re- 
cent mass estimates (see Table 1) have an uncertainty near 2%. The mass estimates 
for Vesta are derived from close flybys of other asteroids with Vesta (Michalak 2000; 
Baer et al. 2008; see home.earthlink.net/~jimbaerl/astmass.txt; Kuzmanoski et al. 2010), the 
effect of Vesta on the Mars ephemeris (Standish 2001; Pitjeva 2005; Konopliv et al. 2006, 
2011; Fienga et al. 2009), the effect of Vesta on the orbit of Eros using the NEAR range data 
(Konopliv et al. 2002, 2011), and a proposed range of values (Pitjeva and Standish 2009). 
The Dawn gravity and mission design studies have adopted a value of GM = 17.5 km 3 /s 2 
for Vesta, picked to be in the middle of the probable range of values (17.3-17.8 km 3 /s 2 ) 
from the latest studies. The recent planetary ephemeris and NEAR estimates of the Vesta 
mass have been improved due to the 10 years of 1-m ranging accuracy to Mars. The NEAR 
Vesta mass estimate (Konopliv et al. 201 1) is an improvement over earlier estimates because 
systematic errors in the Mars perturbation of the Eros orbit are greatly reduced by the several 
orders of magnitude improvement in the Mars ephemeris. The Vesta mass estimates from 
the Dawn mission will be at least 100 x more accurate than the current 2% uncertainty. 

The shape of Ceres was determined to be an oblate spheroid from Hubble images (mean 
radius 476.2 ±1.7 km, Thomas et al. 2005) and from Keck images (mean radius 467.6 ± 
2.2 km; Carry et al. 2008). Although individual volume estimates claim uncertainties near 
1%, the reported values differ by 5%. As with Vesta, the mass of Ceres is also known from 
asteroid close encounters (Michalak 2000; Baer et al. 2008) and from the effect on the Mars 
ephemeris (Standish 2001; Pitjeva 2005; Konopliv et al. 2006, 2011; Fienga et al. 2009) with 
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Table 1 

Selected mass estimates for Vesta and Ceres 


Asteroid 

GM (knrVs 2 ) Mass (1CT 12 M A ) 

Reference 


18.1 ±0.7 

136 ±5 

Michalak (2000) 

17.8 ±0.2 

134.1 ± 1.5 

Standish (2001) 

17.84 ±0.27 

134.4 ±2.0 

Pitjeva (2005) 

18.03 ±0.21 

135.8 ±1.6 

Konopliv et al. (2006) 

17.51 ±0.04 

131.9 ±0.3 

Baer et al. (2008) 

17.9 ±0.4 

135 ±3 

Pitjeva and Standish (2009) 

17.26 ±0.04 

130.0 ±0.1 

Kuzmanoski et al. (2010) 

17.38 ±0.27 

131.0 ± 2.1 

Konopliv et al. (2011), Mars ephemeris 

17.66 ±0.2 

133.1 ± 1.5 

Konopliv et al. (2011), NEAR 

62.4 ±0.2 

470 ±4 

Michalak (2000) 

63.2 ±0.2 

476.2 ± 1.5 

Standish (2001) 

63.08 ±0.27 

475.3 ±2.0 

Pitjeva (2005) 

62.36 ±0.37 

469.9 ±2.8 

Konopliv et al. (2006) 

63.11 ±0.11 

475.5 ±0.8 

Baer et al. (2008) 

62.6 ±0.4 

472 ±3 

Pitjeva and Standish (2009) 

62.10 ±0.43 

467.9 ±3.3 

Konopliv et al. (2011) 


an accuracy of about 2% (see Table 1). The Dawn mission will improve knowledge of the 
mass of Ceres by more than two orders of magnitude in comparison to current estimates. 

The Dawn spacecraft will arrive at Vesta in July 2011 and depart July 2012 (Rayman 
et al. 2011a). During its stay at Vesta, gravity science data will be collected for a total 
of 50 days in two high altitude mapping orbit (HAMO) phases with spacecraft altitudes 
near 700-km and for 70 days in one low altitude mapping orbit (LAMO) science phase 
of the mission with spacecraft altitudes near 200-km. The Dawn spacecraft subsequently 
arrives at Ceres February 2015 and later collects 21 days of gravity data in the higher 
altitude orbit and one month of data during the lower altitude orbit. The mission con- 
cludes at the end of July 2015. The gravity data consists of Deep Space Network (DSN) 
X-band Doppler tracking and optical landmark data. A similar approach of combining 
Doppler tracking and optical landmark data was used for the NEAR mission at Eros. The 
approach will be similar here, except that, instead of using crater centers, the landmarks 
are more generalized control points centered at local topography maps, called maplets that 
cover in a near-uniform manner the visible surface. The topography maps and associated 
landmark network will be constructed using stereo-photoclinometry (SPC) (Gaskell 2005 ; 
Gaskell et al. 2008) and stereo Digital Terrain Maps (DTMs) (Gwinner et al. 2010). 


2 Radiometric and Optical Measurements 

Critically relevant to the Dawn gravity investigation is the telecommunications link between 
the spacecraft and the Deep Space Network (DSN). The right-hand circularly polarized X- 
band link is essentially the mission’s gravity instrument. The uplink signal, for sending 
commands to the spacecraft, is centered at 7179 MHz and the downlink signal, for science 
and engineering telemetry packets is centered at 8435 MHz. Dawn’s navigators and Radio 
Scientists utilize this link in its nominal two-way coherent mode for Doppler tracking in 
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Fig. 1 The angle between the 
Sun and (a) Vesta and (b) Ceres 
as viewed from the Earth for the 
duration of the Dawn mission 
science phases. As the Sun angle 
becomes smaller, the Doppler 
noise increases. The area in red 
indicates the times when the Sun 
angle does not satisfy the 
Doppler accuracy requirement 




order to accurately construct the spacecraft’s motion and small perturbations in its accel- 
erations. In this mode, an uplink signal is derived from the highly stable hydrogen maser 
atomic clock (frequency reference) at the ground station and transmitted to the spacecraft 
which, one-way light time later, locks onto it in the transponder and generates a coherently 
related downlink signal that maintains the phase stability when received back at the ground 
station. The typical resulting noise in the Doppler observables can reach a data quality of 
~0.05 mm/s or better at 10-second integration times. This measure of the gravity science 
data quality depends on several factors that aggregately constitute the experiment’s error 
budget as detailed in Asmar et al. (2005), including the effects of media through which the 
radio beam propagates and the standard calibration techniques for those effects. Two error 
factors, however, are relevant to the planning and execution of the Dawn gravity observa- 
tions: the angular proximity of the radio beam to the Sun and the received signal-to-noise 
ratio (SNR). 

The first factor, the angle between the Sun and Vesta as viewed from the Earth (see 
Fig. 1), bounds the duration of the experiment, since the solar plasma significantly increases 
the Doppler noise at angles under 20°. As an example, Fig. 2 shows how the Doppler ac- 
curacy degrades as Sun angle decreases for the X-band Mars Global Surveyor (MGS) and 
Mars Odyssey missions. The accuracy versus Sun angle is expected to be the same for 
Dawn. Since the Dawn gravity investigation requires the Doppler accuracy for the mission 
to be better than 0.15 mm/s at 10-s samples, the tracking coverage for the gravity experi- 
ment is to occur for Sun angles >20°. The history of the Sun angle for all mission phases in 
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Fig. 2 Doppler noise accuracy for the MGS and Mars Odyssey missions versus Sun angle. Similar accuracy 
is expected for the Dawn mission 


shown in Fig. 1 . The Dawn spacecraft will be first placed into a survey orbit around Vesta for 
17 days beginning August 10, 2011, and the initial Sun angle is near 150°. After a transition 
orbit, the first HAMO phase of 30 days begins September 23, 2011, where the sun angle is 
about 90°. However, the Sun angle decreases to much smaller values for the 70 day LAMO 
phase that starts December 11, 2011. A second 20-day HAMO phase then begins June 8, 
2012 after superior conjunction. With the current mission plan for Ceres, the Sun-Earth- 
Ceres angle is below 20° for the 21 day HAMO orbit from early May to late May 2015, but 
the Sun angle starts at 35° in late June 2015 for the LAMO phase and increases thereafter. 
Note that the actual science phase dates may change due to mission re-plans, but the orbit 
geometry and tracking schedule will remain similar to those that are assumed in this paper. 

The second factor of SNR determines the usability of low gain antennas in the context 
of other instrument data acquisition that can lead to additional gravity data. Although the 
Doppler noise does not decrease with higher SNR, there is a minimum level threshold to 
overcome the phase jitter of the link’s thermal noise. For a carrier frequency / (8435 MHz), 
integration time T (60 s), transponder turn-around ratio G (880/749), and transponder carrier 
loop bandwidth B in Hz, the error contribution a v in units of velocity due to thermal noise 
is, based on Sniffin et al. (2000): 


c / 1 G 2 B 

cry = / h (1) 

2\f2nfT\ Pdown (Pc/No) up 

where pdown is ratio of carrier to noise power in the downlink receiver loop and ( Pc/No) U p 
is the uplink carrier power to noise spectral density ratio, Hz (the latter of which is typically 
not a limiting factor due to the high power transmitters and gain of the ground stations). 

Dawn’s transponder is the Small Deep Space Transponder type, common to most NASA 
deep space missions of the era and proven for excellent performance for Radio Science at the 
~3.6 cm X-band wavelength. The receiver acquires and tracks the uplink carrier via a phase- 
locked loop and produces a voltage controlled oscillator signal whose phase is coherent with 
the uplink carrier. For the coherent downlink transmission mode required for gravity and 
navigation, the exciter derives a coherent X-band carrier from the receiver frequency. When 
coherency is not achievable (e.g., no uplink from a ground station), the downlink carrier is 
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derived from the auxiliary oscillator and the Doppler noise of this internal frequency source 
is too high to be useful to gravity science. 

The telecommunications subsystem is equipped with a traveling-wave tube amplifier 
(TWTA) that raises the transponder’s output signal level to 100 Watts prior to radiating 
through the antenna in order to meet the telemetry highest bit rate requirement at the furthest 
distances of the mission. The High Gain Antenna (HGA) is a 1.52-m diameter paraboloid 
with a nominal downlink bore sight gain of 39.6 dBi (Taylor 2009). It is fixed mounted to 
the +X panel of the spacecraft, thus requiring spacecraft pointing to Earth. There are also 
three Low Gain Antennas (LGA) whose peak gain is ~6 dBi, mounted along the +x, +z, 
and — z axes, that sacrifice gain for beam width in order to provide a uniform coverage over 
various spacecraft orientations for purposes of low data rate communications, especially 
during engine thrusting, and during safe mode. It is expected that the spacecraft’s HGA will 
be used for the planned gravity observations, providing sufficient SNR in the link to meet 
the Doppler noise requirements. The LGAs can also be used to provide additional data when 
the HGA is not pointed to Earth to satisfy the needs of other instruments. To achieve the re- 
quired minimum 0.15 mm/s at 10 s integration, the received carrier to noise power must be 
8.6 dB-Hz and 22.5 dB-Hz to reach the goal of 0.03 mm/s. Analysis of Vesta LAMO orbit 
geometry showed that 16 dB-Hz are expected from an LGA to a 34 meter diameter DSN 
station. An LGA link to a 70-m station provides an additional 6-8 dB thus reaching the goal 
with sufficient margin. 

The nominal amount of tracking for Dawn during the gravity phases of the mission is 
three 8-h (first 6 weeks) or three 5.5-h (next 4 weeks) passes using the HGA per week, 
one LGA pass per day using the 70-m DSN antennas, and additional continuous tracking 
using the LGAs with the DSN 34-m antennas. Expected results are given below for two 
scenarios — one 5.5-h pass per day and continuous tracking. Actual results are expected to 
be between the two. The ground track spacing for the HGA plus 70-m LGA tracking satisfies 
the 30-km and iOO-km spacing requirement at Vesta and Ceres, respectively. In addition to 
Doppler tracking, one to two meter accuracy range data will be acquired, and this is essential 
to eliminating Vesta or Ceres ephemeris errors from the gravity solution. 

Science images will be acquired and used for optical navigation at all mission phases 
(Polansky et al. 2011). The main optical observables will be landmark vectors from the 
center of the body to the surface. Landmarks are control points at the center of a local digital 
elevation model or maplet. A maplet is defined by its control point, its local coordinate 
system and a local grid, typically 99 x 99 pixels. Each grid-point contains the local elevation 
relative to the control point as well as a relative albedo value. 

The details of landmark-based optical navigation for Dawn have been presented else- 
where (Mastrodemos et al. 2011). Here we give a brief summary. Each maplet and its asso- 
ciated landmark vector is constructed as follows: (1) combine images at varying view angles 
and incidence angles, with a minimum of 3 images needed, (2) from each image set extract a 
template of a predefined grid scale and size and normally-project it onto the a-priori asteroid 
surface (a-priori shape model), (3) for each image template construct a brightness model, a 
step which includes an assumption about the reflectance model of the asteroid, (4) perform 
a least-squares fit to the brightness model and estimate at each maplet pixel the local slopes 
and relative albedo, (5) integrate the slopes to construct the heights. 

Once a number of maplets have been constructed a global geometric solution is per- 
formed (see Lig. 3), that solves for the landmark vector (one per maplet), the camera point- 
ing and the spacecraft-asteroid center relative vector, the latter in a kinematic sense, i.e., on 
a per-image basis. Subsequently the existing maplets are re-derived, this time starting with 
better initial conditions of camera pointing and landmark vectors. Additional data types that 
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Fig. 3 A simulated Vesta Survey 
image of the south pole region. 
The white square circles mark the 
locations of the landmark vectors 
that are included in this image 



contribute to this solution are the maplet-to-maplet overlap vectors and the limb apparition 
of maplets. This process of maplet construction, global geometric solution and outlier rejec- 
tion iterates at least 2, 3 times with each data set, prior to adding new images and/or new 
maplets. 

At every mission phase a network of landmarks that correspond to tightly overlap- 
ping maplets is constructed. The landmark network begins with relatively coarse resolu- 
tion maplets and sparse spacing of landmark vectors from the approach observations and 
it is subsequently refined at the next mission phase. The optical measurements consist of 
the image pixel, line location of each landmark and of a Cartesian vector from the center 
of mass to the surface (landmark vector). Assuming that a number of images, typically 5 
or 6, at varying incidence and emission angles are included in each landmark, to first order, 
both the accuracy of the landmark pixel, line location in the image and the accuracy of the 
landmark vector will depend on the image resolution and the maplet scale. Typically, at the 
beginning of each processing phase, the maplet scale is coarser by up to twice that of the 
image scale and at the end of each processing phase it is comparable to that of the image 
resolution. If a large number of images, larger than ~10, has contributed for the construc- 
tion of each maplet, then the problem of the local topography estimation can become highly 
over-determined. In this case, as with the Survey observations, maplets up to 1/2 the image 
scale can be constmcted, allowing even higher landmark accuracy, although in this case the 
scaling may not be linear with maplet resolution. 

Typical 1 — a accuracy estimates at the end of each processing phase are for a post-fit 
image pixel, line residuals of ~ 1/3 image pixels and for the landmark vector ~ 1/3 of the 
corresponding maplet scale (km/pixel) per Cartesian component. For example, at the last 
rotational characterization at a range of ~5000 km we expect pixel, line post-fit residuals of 
~ 150 m and landmark vector accuracy of ~ 1 25- 1 50 m. According to the current mission 
plan at Vesta the approximate number of images, both for optical navigation and science, and 
the respective number of landmarks and optical data errors for each mission phase are given 
in Table 2. Note that the LAMO images are acquired all at a nadir direction and although 
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Table 2 Expected number of images and performance for optical landmark tracking for each Dawn mission 
phase at Vesta 


Mission 

Phase 

Total # of images 
(mean resolution) 

Approximate 
total # of 
landmarks 

Average # of 
landmarks 
per image 

Landmark 
vector 1-er 
Position error 

Landmark pixel, 
line post-fit 
residual in meters 

Survey 

520 (280 m) 

12000 at 
150 m/pixel 

>600 

140 m 

140 m 

HAMOl 

2400 (60 m) 

120000 at 
50 m/pixel 

50 

40 m 

40 m 

HAM01+ 

HAM02 

3900 (60 m) 

120000 at 
50 m/pixel 

30 

40 m 

40 m 

LAMO 

2000 (25 m) 

120000 at 
50 m/pixel 

30 

40 m 

40 m 


they allow an almost complete coverage of the asteroid surface, they are regarded more as 
supplementary to HAMO images, rather than a complete data set in its own right. 


3 Spacecraft Models 

The major features of the Dawn spacecraft including the spacecraft coordinate frame, sci- 
ence instruments, and antenna locations are shown in Fig. 4. The high aspect-ratio solar 
panels are a distinctive feature of the Dawn spacecraft and represent a larger surface area 
(36.5 m 2 , ~50% increase in area to mass ratio) than other spacecraft such as MRO, MGS, 
Odyssey and NEAR because of the power required for the ion engines. The central ion en- 
gine is located along the — z thrust direction, and the other two ion engines are located in 
the x-z plane and are canted about ±40° from the — z direction. No gravity data is pro- 
cessed when the ion engines are in operation. The solar panels, which can rotate along the 
y-axis, always points near normal to the Sun direction during all phases of the mission. 
The spacecraft attitude is determined by the power- steering algorithm, which computes the 
attitude given a spacecraft body-fixed aim vector (e.g. instrument deck or +z direction for 
science observations and +x direction for HGA tracking) and a target vector (e.g., Vesta 
for science observations and Earth for HGA tracking). The body-fixed aim vector is first 
aligned with the nadir or Earth pointed target vector, during most of science phases, and the 
spacecraft rotates about the aim vector axis and finds a solution which in order of priority 
puts first the -x-axis (Visible InfraRed Spectrometer or VIR, radiators, and star trackers) 
or second the +z-axes (Gamma Ray and Neutron Detector or GRaND, Framing Camera, 
and gyros) away from the Sun. If multiple solutions exist, the solution with the least angular 
displacement is selected. The spacecraft attitude control subsystem provides attitude deter- 
mination and control during all mission phases and consists of the following elements: two 
star trackers looking in the — x direction, an inertial reference unit and four reaction wheel 
assemblies. The latter are required to maintain stability in the spacecraft motion to mini- 
mize non-gravitation acceleration, which would be another noise factor for gravity science. 
The quaternions for the spacecraft attitude and solar array orientation with respect to the 
spacecraft frame will be obtained from the Planetary Data System (PDS) archive (Semenov, 
private communication, 2010). 

During gravity science phases, a body-fixed high-gain antenna (HGA) will be used for 
tracking when it’s available; otherwise, low-gain antennas (LGAs) will be used. When the 
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Fig. 4 The Dawn spacecraft 
with location of science 
instruments and antennas. The 
spacecraft coordinate frame is 
shown; its origin is at the 
baseplate of the central ion 
engine 



Flight configuration 


Solar arrays 
(rotates about y) 


HGA is in use, the HGA will be pointing toward the Earth, and phase center location with 
respect to the spacecraft center-of-mass remains constant in the Earth direction. Thus, the 
spacecraft attitude is not needed in the HGA Doppler measurement model. However, when 
LGAs are used, the spacecraft will be slowly changing its attitude with respect to Earth to 
maintain nadir point, which requires the LGA motion with respect to the center of mass in 
the Doppler measurement model. Note that there are three LGAs as shown in Fig. 4. 

Figure 5 shows the current reference trajectory designs for HAMO and LAMO at Vesta 
and Ceres. The center bodies are drawn to scale to reflect the orbit sizes (based on cur- 
rent shape estimates). The major forces acting on the spacecraft are the accelerations due to 
gravity and solar radiation pressure. Figures 6 and 7 show the expected accelerations during 
HAMO and LAMO mission phases. Note that CBP refers to the central body point- mass 
acceleration, OBL refers to the higher-order gravity acceleration (non-point mass accelera- 
tion), and SRP refers to the acceleration due to solar radiation pressure. The trajectory was 
integrated using JPL’s DE421 ephemerides for positions of the planets. 

Since Vesta is a differentiated body (e.g., Keil 2002), a simple three-layered model was 
used as the nominal field for generation of the reference orbit (similar to models tested by 
Tricarico and Sykes 2010). The 15th spherical harmonic degree gravity field model consists 
of a spherical 150-km radial core of density 5.8 gm/cm 3 at the center-of-figure, an ellip- 
soidal mantle (250 x 250 x 227-km and 3.0 gm/cm 3 ) and a crust (2.7 gm/cm 3 ) from the 
shape model (Thomas et al. 1997) with variable thickness (zero near the southern crater 
and ~20-40-km elsewhere). Multiple models were tested and the accuracy of the gravity 
field that will be determined at Vesta is not that sensitive to the nominal field chosen for 
the simulation (various core densities and location or even a uniform density Vesta). The 
gravity acceleration for Vesta in Fig. 6 uses the nominal 3-layer model. The accelerations 
from an assumed uniform density gravity field model are about 20% higher and represent 
a probable upper bound of all various internal models. The Ceres gravity acceleration used 
a 20th degree gravity field which was sampled based on a Kaula rule, i.e., K Cer es = 0.0086 
(see below). 

The solar pressure acceleration is computed using the appropriate specular and diffuse 
reflectivity coefficients for a multiple flat-plate box-wing spacecraft model. Although the 
acceleration due to solar radiation pressure is several orders of magnitude smaller than the 
gravity acceleration (and 50% higher than typical due to a larger solar panel), it is one of the 
limiting error sources for high-precision determination of the gravity field and pole location. 
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-HAMO 

-LAMO 


Sun 

Direction 



Pole 

Vector 


Fig. 5 Views of the HAMO and LAMO orbits at (a) Vesta as viewed from the ecliptic north pole and Sun 
and (b) Ceres as seen from the ecliptic north pole. The directions of the Sun and Earth are also shown 


Another important force that will affect the gravity science is the angular momentum 
dump (AMD), which stabilizes the on-board wheel speeds using hydrazine thrusters. The 
largest expected magnitude of AMDs during HAMO and LAMO phases is near 10 mm/s 
with a duration of 10-20 minutes. Also note that AMD magnitude will be higher during 
LAMO than HAMO since the gravity gradient force increases as the orbit radius decreases. 
The magnitudes are not as important as the AMD frequency since corrections to the AMD 
amplitudes will be estimated as a delta-V maneuver. The longer the duration between the 
AMDs, the better the gravity information in the data arc. The current science performance 
presented below assumes an AMD maneuver every two days for both HAMO and LAMO 
(close to the planned three times per week for LAMO but more conservative for the expected 
once per week in HAMO). During the LAMO phases, the current nominal plan assumes that 
there will be an orbit maintenance maneuver (OMM) every 7 days. 


4 Gravity Determination Technique 

The gravity fields of Vesta and Ceres will be determined through the gravitational pertur- 
bations of the orbit of the spacecraft. The technique and software (Orbit Determination 
Program or ODP, Moyer 1971) are the same as that used for previous planetary and small 
body gravity analysis (e.g., Konopliv et al. 2011 for Mars and Konopliv et al. 2002 for 
Eros). The equations of motion and variational partials for the estimated parameters are nu- 
merically integrated (variable order Adams method, Krogh 1973) for the nominal motion 
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Fig. 6 The accelerations on the Dawn spacecraft at Vesta due to the central two-body point mass (CBP), 
the gravitational oblateness (OBL) and solar radiation pressure (SRP). The constant density Vesta gravity 
model was used to compute the oblateness acceleration. The accelerations are given the radial ( R ), along the 
velocity (T) and normal to the orbit plane ( N ) directions 


of the spacecraft in the International Celestial Reference Frame (ICRF, within <100 mas of 
Earth Mean Equator of 2000; Folkner et al. 1994). All forces acting on the spacecraft are 
accounted for including gravity, solar radiation pressure, thrusting (mostly for despin of the 
angular momentum wheels), and point mass contributions from the Sun and planets. The 
DSN Doppler and landmark observations are then processed (Moyer 2000) to obtain the 
linearized observation equations (e.g., Tapley et al. 2004). 
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Fig. 6 ( Continued ) 


The spacecraft state and other parameters are estimated using a square root information 
weighted least squares filter (or SRIF, see Lawson and Flanson 1995, and Bierman 1977). 
The observations are divided into time spans that we call data arcs. For the Dawn mission we 
expect arc length to be about 4 days for LAMO, but longer arc lengths will be investigated 
to the point where gravity results are not degraded by modeled non- gravitational forces. The 
estimated parameters consist of local parameters specific to the data arc such as spacecraft 
position and velocity, velocity corrections for angular momentum wheel despins, and solar 
pressure coefficients and global parameters common to all data arcs such as gravity, pole 
location and rotation rate. The global parameters are determined from merging the global 
SRIF data array from each arc similar in concept to the partitioned normal matrix method 
(Kaula 1966). 

The gravitational potentials of Vesta and Ceres are modeled by a spherical harmonic 
expansion in the body-fixed reference frame with normalized coefficients (C nm , S nm ) that is 
given by (e.g., Heiskanen and Moritz 1967; Kaula 1966) 


GM GM / R e . 

U = — H — 2_^ 2^1 T ) ^imCsin 4>i at )\C nm cos (mk) + S nm sin (mX)], (2) 

n= 1 m= 0 ' 


where GM is the gravitational constant times the mass of the central body, n is the degree, 
m is the order, P nm are the fully normalized associated Legendre polynomials, R e is the 
reference radius of the body (265-km for Vesta, 467-km for Ceres), 4>i at is the latitude, 
and A is the longitude. The gravity coefficients are normalized such that the integral of 
the harmonic squared equals the area of a unit sphere, and are related to the unnormalized 
coefficients as (Kaula 1966; Lambeck 1988) 



in - m)\(2n + 1)(2 - 8 0m ) 


-|!/2 


in + m)\ 



( 3 ) 
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Fig. 7 The accelerations on the Dawn spacecraft at Ceres due to the central two-body point mass (CBP), 
the gravitational oblateness (OBL) and solar radiation pressure (SRP). A sample 20th degree harmonic field 
based upon the expected Kaula power rule was used for the oblateness acceleration. The accelerations are 
given the radial (R), along the velocity (T) and normal to the orbit plane (N) directions 


The degree one coefficients are zero, since the origin of the coordinate system is chosen to 
be the center of mass. The zonal coefficients are given by J„ = — C„o- 

The RMS magnitude spectrum of the gravity field tends to follow a power law of the 
form (Kaula 1966) 


Springer 


473 


Reprinted from the journal 


A.S. Konopliv et al. 



Fig. 7 ( Continued ) 




t:=o(cl + S 2 nm ) ~ K 


2n + 1 


( 4 ) 


where the constant K is ~10 -5 for the Earth (values for Vesta and Ceres are discussed 
below). A constraint of this form is often applied to the gravity coefficients in order to limit 
excessive power in the poorly determined coefficients of higher degree. The coefficients are 
biased to zero with an uncertainty given by the power law. To avoid suppressing gravity 
information in the better- known coefficients, the constraint is sometimes only applied to 
parts of the spectrum where M„ exceeds the power law. 

Combining spherical harmonic coefficients from the anticipated shape models for Vesta 
and Ceres (Raymond et al. 201 lb) will allow long- wavelength estimates of interior structure 
(cf. McCord et al. 2011; Zuber et al. 2011) posed in terms of crustal and elastic thicknesses 
or mantle density variations (Neumann et al. 1996; Wieczorek 2007). 


5 Moments of Inertia from Gravity, Wobble and Obliquity 

In addition to determining the mean densities of Ceres and Vesta, the Dawn mission may be 
able to obtain estimates of their moments of inertia. Together, the mean density and mean 
moment of inertia provide two integral constraints on the radial density profile, and allow 
better determination of the internal density structure (Bills and Rubincam 1995). 

There is a close relationship between the inertia tensor of a body and the degree two 
coefficients in a spherical harmonic expansion of the gravitational potential (Soler 1984). 
However, there are only 5 degree two coefficients and 6 independent terms in the inertia 
tensor. Thus, even if we have a highly accurate gravity field, we will still lack sufficient 
information to directly determine the moments of inertia. We now discuss three possible 
means of estimating the remaining constraint. 
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5 . 1 Moment from gravity 


If the coordinate axes are chosen to coincide with principal inertial axes, the inertia tensor 
diagonalizes, three of the potential coefficients vanish and the remaining two unnormalized 
gravity coefficients, which we denote J 2 and C 2 . 2 , have a simple relationship to the principal 
moments A < B < C 


h 

MR 1 = 

C-(A + B)/2 

_C 2 ,2_ 

e 

( B - A) /A 


(5) 


where M and R are the mass and mean radius of the body, respectively. 

If a rotating body is in hydrostatic equilibrium, there is a unique relationship between the 
imposed rotational potential and the induced gravitational potential. The mean moment of 
inertia 


I = (A + B + C)/3 


( 6 ) 


is not changed by tidal or rotational deformations (Rochester and Smylie 1974), but the 
principal moments are individually changed. In much of what follows, it will be convenient 
to use dimensionless versions of the principal and mean moments of inertia, which are given 
by 
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~A~ 
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B 
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C 
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I 



We can write the rotational perturbation as (Bills 1995) 


(7) 
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-1 

B 

= 1 + 
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kor /?/ 
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(8) 


where 00 is the rotation rate, G is the gravitational constant, and k is a secular Love number 
(Munk and MacDonald 1960). The corresponding values of the degree two gravitational 
potential coefficients are 


r h 1 


1 

_C 2j2 _ 


0 


2d3 


kco'R L 
3 GM 


(9) 


Measurement of Ji will provide an estimate of the secular Love number, and the ratio 
of C 2 , 2 // 2 will provide an estimate of the proximity to hydrostatic equilibrium. If the hy- 
drostatic assumption is justified, the secular Love number can be used to estimate the polar 
moment of inertia from the Radau-Darwin relation (Radau 1885; Darwin 1899) 


C 2/ 2 [ 4=k \ 

MR] 3 \ 5V 1 +k) 


( 10 ) 


This relation holds exactly only for uniform density fluids, but is a reasonable approximation 
for radially stratified bodies (Nakiboglu 1982). 

Current estimates of the moments of inertia of the Galilean satellites of Jupiter all depend 
upon this approximation (Hubbard and Anderson 1978; Anderson et al. 1996, 1998a, 1998b, 
2001). There are two primary differences between those cases and the Dawn encounters with 
Ceres and Vesta. The Galileo estimates are all derived from multiple fly-by opportunities, 
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rather than going into orbit about the bodies of interest. Also, since the Galilean satellites 
are all in synchronous rotation, they have appreciable tidal contributions, which add to the 
J 2 values, and make the expected € 2,2 values non-zero. 

5.2 Moment from wobble 

If the spin pole of either Ceres or Vesta exhibits a free wobble, or misalignment with the 
axis of greatest inertia, that would be informative in two ways. First, it would suggest a very 
recent excitation event, and second, it would provide a means of estimating the moment of 
inertia. 

The rotational configuration of lowest energy, for a given angular momentum, consists 
of steady rotation about the axis of greatest moment of inertia. If the body is rotating about 
another axis, it will exhibit a wobble. That is, the instantaneous rotation vector does not 
maintain a fixed orientation relative to the body. The free wobble angular rate w is, like the 
spin pole precession rate, also determined by the spin rate and the moments of inertia. It is 
given by (Smith and Dahlen 1981) 


Bills and Nimmo (2011) have estimated wobble periods of 5.5 and 1.2 days for Ceres and 
Vesta, respectively. If the wobble period could be accurately determined, it would provide 
one additional constraint, and thereby allow the moments of inertia to be estimated. 

There have been numerous analyses of wobble damping, as applied to asteroids and 
satellites (Bums and Safronov 1973; Harris 1994; Black et al. 1999; Sharma et al. 2005; 
Samarasinha 2008). The basic idea is that the body is deformed, in response to the rotational 
potential, and that part of this rotational bulge follows the instantaneous equator as it moves 
over the surface of the body. The continuing deformation dissipates energy and drives the 
body toward the lowest energy state of principal axis rotation. Bills and Nimmo (2011) 
estimate wobble damping times for Ceres and Vesta of 440 and 120 years, respectively. 
Unless these estimates were based on assumptions that are overly simplistic, or there has 
been a very recent excitation event, a wobble would not be expected to be observed. 

5.3 Moment from obliquity 

If, as expected, the spin poles of Ceres and Vesta are fully damped, then measurements of 
their obliquities will also provide information about their moments of inertia. It would also 
provide information about the amount of dissipation within these bodies. This approach has 
recently been applied to Titan (Bills and Nimmo 2008) and more recently Ceres and Vesta 
(Bills and Nimmo 2011). 

The moments of inertia are most directly linked to rotational dynamics and the standard 
astronomical estimates of moments of inertia arise from determination of spin pole preces- 
sion rates. That is how the moments of inertia of Earth (Williams 1994) and Mars (Folkner 
et al. 1997) are known. However, the expected precession rates for Ceres and Vesta are small 
(few arcsec/year) and there is little hope of directly determining precession during the short 
visits by Dawn. 

An alternative exists, which is that, if the spin pole is fully damped, then its instantaneous 
orientation, relative to the orbit pole, will be diagnostic of the precession rate. The Moon 
is a classic example of such a body. The orbit pole precesses about the ecliptic pole, and 
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the spin pole precesses about the orbit pole. However, due to the achievement of a steady, 
damped configuration, the spin pole, orbit pole, and ecliptic pole remain coplanar. Such a 
configuration is now known as a Cassini state (Colombo 1966; Peale 1969; Henrard and 
Murigande 1987). 

The condition for this co-planar precession, in steadily precessing nearly circular orbits, 
can be written as (Colombo 1966) 


where i is the inclination of the orbit pole to the invariable pole, s is the obliquity or sepa- 
ration of spin and orbit pole. The parameter /l is the rate of orbit plane precession, and a is 
a spin pole precession rate. It is determined by the moments of inertia of the body, and the 
relative rates of orbital motion n and spin s. It has explicit form 


When the obliquity s and gravity coefficients and C 2 .? are all known, the constraint 
equation (12) can be solved for the polar moment of inertia, without any assumption of 
hydrostatic equilibrium. However, it does depend upon an assumption that the spin pole is 
damped. 

How likely is it that the spin poles of Ceres and Vesta are damped? For small asteroids, 
non- gravitational effects can have appreciable influence on rotation rate and direction. For 
example, many members of the Koronis family of asteroids have non-random spin rates and 
directions (Slivan 2002; Slivan et al. 2003, 2008), counter to the expectation for purely colli- 
sional evolution. This quasi-alignment has been attributed to thermal torques (Vokrouhlicky 
et al. 2003, 2006; Bottke et al. 2006; Scheeres and Mirrahimi 2008). 

For larger and more nearly spherical bodies like Ceres and Vesta, the mechanisms of 
excitation and damping are less clear. If the near-surface of Ceres is ice-rich (Thomas et 
al. 2005), its viscoelastic relaxation (Maxwell) time will be short, perhaps even compara- 
ble to the rotation period, resulting in enhanced dissipation. Peale (2005) has argued that, 
for Mercury, where tidal damping is reasonably effective, the dominant obliquity damping 
mechanism is likely viscous core-mantle coupling. If large asteroids have (or had) fluid re- 
gions within them, their spin poles could be fully damped. It has been argued that at least 
one large asteroid, the Angrite parent body, possessed a fluid core early in its history (Weiss 
et al. 2008). 

One means of testing for an un-damped obliquity would be to look for a free wobble, 
or non-principal axis rotation (Munk and MacDonald 1960; Lambeck 1980; Chen et al. 
2009). The wobble period is short compared to the spin pole precession period, and thus 
the wobble damping time is short compared to the obliquity damping time. A hypothetical 
excitation event, such as a relatively recent impact (Peale 1976), could have produced both a 
free wobble and a free obliquity. If no wobble is found, it would be permissive of a damped 
obliquity. However, since the obliquity damping time is expected to be much longer, absence 
of a wobble is no guarantee of a fully damped obliquity. If, on the other hand, an un-damped 
wobble is found, it would strongly suggest recent excitation. While such a wobble would 
presumably indicate that the obliquity is not fully damped, measuring the wobble period 
would provide an alternative way of determining the moment of inertia. 

Another means of testing the hypothesis of a damped spin pole is to examine the az- 
imuthal orientation of the spin pole, relative to the orbit pole. As is shown in Bills and 
Nimmo (2011), if the spin pole is fully damped, its present orientation is fully specified by 
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a single parameter, which is the spin pole precession rate. If the observed spin pole does not 
lie on the curve specified by this model, that would imply there is still a finite free term in 
the spin pole motion, and the spin pole is not fully damped. 


6 Expected Gravity Science Results at Vesta 

Although there is no direct a priori information for the Vesta gravity field, the expected 
plausible range of possible gravity fields (e.g., Tricarico and Sykes 2010) can be inferred 
from the Vesta shape model as determined by the Hubble images (Thomas et al. 1997) and 
the scaling of the gravity coefficients to an expected power law K Ves ta/n 2 (Kaula 1966). 
Figure 8 shows the gravity field without J 2 on the surface of a sphere (290-km or the maxi- 
mum equatorial radius to assure convergence near the Brillouin sphere) assuming our 3-layer 
model for Vesta. The range in surface acceleration is ~2450 mGals. The equivalent range 
from a uniform density model is ~30% greater. 

Lack of convergence within the Brillouin sphere (the smallest sphere that encloses Vesta) 
is a limitation of the spherical harmonic formulation. A map of the gravity near the surface 
of Vesta requires alternate representations of the gravity potential such as ellipsoidal har- 
monics (Garmier and Barriot 2001, and Garmier et al. 2002) that we will investigate as well. 
The convergence of the spherical harmonics is tested in similar fashion to Garmier and Bar- 
riot (2001). The acceleration on an ellipsoid near the Vesta surface was computed using a 
uniform density polyhedral gravity model (Wemer and Scheeres 1997) and the equivalent 
spherical harmonic representation. For a 20th degree harmonic field and radii of 265-km 
or greater (or between latitudes of ±50° or ~75% of the surface), the RMS difference was 
13 mgals (<1% of Fig. 8) and a maximum of ~ 170 mgals (or 8%). However, gravity differ- 
ences were substantially greater (~700 mgals) near the poles with a radius of 243-km and 
require the ellipsoidal harmonic representation. 
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Fig. 8 The accelerations on a 290-km sphere for the 3-layered Vesta gravity model, a 15th degree spherical 
harmonic gravity model based upon the shape of Thomas et al. (1997). The acceleration due to J 2 is removed 
and the maximum and minimum accelerations are 1424 and —1024 rnGals, respectively 
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Fig. 9 The gravity RMS magnitude spectrum for Vesta simulations of one-month duration. The following 
spectrums are shown: the expected Kaula power law, gravity from shape assuming a constant density, gravity 
from shape with a 3-layer model, expected errors from the HAMO science phase using DSN tracking only, 
expected errors from the LAMO science phase using DSN data only, optical landmark tracking only, and DSN 
data together with optical landmark tracking. The amount of tracking included in the DSN data is 5.5 h/day 


The expected power in the spherical harmonic coefficients are given by scaling the 
Earth’s power law by (Kaula 1963) 


K 


Vesta 


= K 


Earth 


( M Earth \ ( R Vesta 


\Mv e sta ) \R 


'Earth 


(14) 


which gives K Vesta = 0.011 for GM Vesta = 17.66 km 3 /s 2 (Konopliv et al. 2011) and 
R Vesta = 265 km (Thomas et al. 1997). For Vesta, the gravity from shape contributions 
at low degree (n <~ 10) are greater than the expected power law, but drops off at higher de- 
gree probably due to loss of image resolution (Fig. 9). The gravity spectrums from various 
layered and shifted core models are similar in size (Tricarico and Sykes 2010). 

After the Vesta Survey orbit at a radius of 2700-km, Dawn will be placed in HAMO- 
1 for 29 days for the first acquisition of gravity data during September 2011. The orbit 
is polar and circular with a radius of 950-km and an orbit period of about 12 h. At this 
altitude, the GM of Vesta is determined to ±0.02 km 3 /s 2 (lOx formal uncertainties), and 
the low degree harmonics of the gravity field are significantly determined through harmonic 
degree 6 (Fig. 9). These HAMO results are from DSN data only with assumed 5.5 h/day of 
X-band 0. 1 mm/s 60-s samples. 

Next, after an orbit transition period, the Dawn spacecraft begins FAMO at Vesta on De- 
cember 12, 2011 in a polar circular orbit. With an orbit semi-major axis of 460-km and orbit 
period of 4 hours, the average altitude of the Dawn spacecraft over the Vesta surface (mean 
265-km radius, Thomas et al. 1997) is about 200 km. Accounting for the shape of Vesta, the 
average altitude over the north and south pole is about 200-km and 250-km, respectively. 
The altitude difference is due to the shape of Vesta (25-km and 45-km below the refer- 
ence radius for the north and south pole, respectively) and orbit eccentricity (30 km higher 
in the south). For a given latitude, the altitude of the spacecraft can vary within a 60-km 
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range due to the strong irregular gravity field of Vesta, which is currently derived from the 
3-layer Vesta gravity model. Tricarico and Sykes (2010) have shown the altitude variations 
are similar for other gravity models including constant density from shape (Thomas et al. 
1997). 

The expected quality of the LAMO gravity field determination at Vesta was investigated 
using 28-days of DSN Doppler tracking and optical landmark images prior to March 1, 
2012. The LAMO radiometric (DSN) only, optical only and radiometric plus optical land- 
mark tracking gravity error spectrums are displayed in Fig. 9. The gravity performance 
results assume 5.5 h/day X-band tracking with 0.1 mm/s accuracy at 60-s samples (easily 
achieved with planned HGA and LGA with DSN 70-m antennas), and the optical data as- 
sumes 20 images over 20 minutes once per day with 1 .0 to 0.4 pixel resolution for an average 
of 30 tracking landmarks per picture. With both DSN and optical data and orbit variability 
assumptions, the spherical harmonic gravity field can be determined to about degree and or- 
der 20 (or 41 -km half wavelength resolution at the surface of Vesta). These results provide 
ample margin for the 90-km or spherical harmonic degree 9 gravity requirement. In addition, 
the GM determination is at least 10 x better than the HAMO expected results. Uncertain- 
ties for periodic signatures in C 21 and Vi (10 -7 , and consistent with the DSN plus optical 
solution of Fig. 9) indicate that wobble detection is possible for signatures on the order of 
10-cm. The gravity information is primarily from the DSN Doppler data. The optical data, 
however, noticeably helps the long wavelength gravity field and orbit determination, as was 
the case for the NEAR mission (Konopliv et al. 2002). 

The Dawn mission at Vesta will also acquire continuous DSN tracking during LAMO 
using the LGA on the spacecraft to the 34-m DSN stations. Figure 10 demonstrates the 
improvement for this scenario with 28 days of tracking assuming the data quality from 
the low gain is comparable to the high gain antenna. The expected resolution improves to 
harmonic degree 27 (or 30-km on the surface). A surprising RMS uncertainty curve results 
from this case, and is probably tied to large variations in the orbit. Plots of gravity coefficient 
uncertainty versus order reveals better determined higher order terms (near sectoral) for 
coefficients of low degree (e.g., n = 5). For higher degree terms, the low order (near zonal) 
coefficients are better determined. Expected results should be between the two cases of 
5. 5 -h/day and continuous tracking. 

The orientation of the Vesta pole and its rotation rate are also determined as part of the 
solution for the gravity field through the dynamics of the spacecraft orbit. This solution is 
independent from the imaging team determination (Raymond et al. 2011b). The uncertain- 
ties in right ascension and declination of the pole (Seidelmann et al. 2007) are determined to 
better than 0.01° with DSN data only, whereas inclusion of the landmark tracking improves 
pole location by at least 5 x (see Table 3). In contrast, the landmark tracking significantly im- 
proves (1000 x ) the rotation rate solution versus DSN only solutions and improves ( 10 x ) the 
current knowledge (Drummond et al. 1988). As for previous gravity studies (e.g., Konopliv 
et al. 2006), the quoted uncertainties are 10 x formal uncertainties to give realistic uncertain- 
ties that account for the red noise characteristics (power spectrum /~ 8/ 3 ) of the Doppler data 
(Woo 1975) and other possible systematic errors such as modeled non-gravitational accel- 
erations on the spacecraft. Any possible free wobble can be detected through time varying 
changes in the C 21 and V| gravity coefficients. The size of a detectable wobble for Vesta 
with period of 1.2 days (see above) is on the order of 20-cm on the surface (r = 265 km; see 
Table 3). 
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Fig. 10 The Vesta gravity error spectrums for the case with continuous (24 h) DSN tracking versus (a) har- 
monic degree and ( b) harmonic order. Continuous tracking is an option for Dawn using the low gain antenna. 
This option improves the resolution from harmonic degree 19 to degree 27, and reduces the uncertainty of 
the determined coefficients by about an order of magnitude. The errors versus harmonic order are shown 
individually for degrees 5, 10, 15, and 20 


7 Expected Gravity Science Results at Ceres 

The gravity field at Ceres is expected to be much smoother than Vesta because of small 
excursions (<10-20 km) from its nearly hydrostatic shape (Thomas et al. 1997; Carry et 
al. 2008) and lower density. In addition to the GM, the only a priori estimates for the Ceres 
gravity field come from the scaling of the Kaula power law ( K Cer es = 0.0086, using Rceres = 
476 km, GMceres = 63.2 km 3 /s 2 ) and predictions of the Ji value coming from the measured 
flattening and assumptions of hydrostatic equilibrium and internal structure. The values in 
unnormalized J 2 range from 0.016 for more differentiated models with a core to 0.032 for a 
homogeneous body (McCord and Sotin 2005; Zolotov 2009). 

The Ceres survey orbit with a radius of 6400 km begins February 18, 2015 and lasts for 

8 days (1.75 orbits). After a transfer phase, the HAMO mission phase begins May 2, 2015 
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Table 3 Vesta and Ceres expected uncertainties in orientation from the Dawn mission. The uncertainties 
are for only the LAMO mission phase and the formal uncertainties are scaled by lOx to give more realistic 
values. The current values used as given as well as a recent pole solution for Vesta 

Asteroid 

R. A., a 

(deg) 

Dec., S 
(deg) 

Rot. rate, u> 
(deg/day) 

Wobble 
amp. (cm) 

Notes 

Vesta 

±0.0040 

±0.0057 

±0.00043 

±330 

DSN 5.5 h/day 


±0.0007 

±0.0008 

±0.0000003 

±300 

DSN 5.5 h/day ± 
Optical Landmarks 


±0.0002 

±0.0003 

±0.0000001 

±24 

DSN 24 h/day ± 
Optical Landmarks 


301 ±5 

41 ±5 

1617.332776 ± 
0.000001 


Current values 
Thomas et al. (1997) 
Drummond et al. (1988) 


305.8 ±3.1 

41.4 ±1.5 



Recent Values 
Li et al. (2011) 

Ceres 

±0.098 

±0.0071 

±0.00075 


DSN 5.5 h/day 


291 ±5 

59 ±5 

952.15320 ± 
0.00021 


Current values 
Thomas et al. (2005) 
Chamberlain et al. (2007) 



Harmonic Degree 

Fig. 11 The gravity RMS magnitude spectrum for the Ceres one-month simulation. The following spectra 
are shown: the expected Kaula power law, the HAMO science phase expected results, and the LAMO science 
phase expected results. These results are for DSN tracking only at 5.5 h/day. Optical landmark tracking is not 
part of the simulation at this time, but it is expected to improve results in similar fashion as with Vesta 


and lasts for 21 days in an orbit with radius 1800-km and a period of 16.8 hours. Since 
the HAMO mission phase occurs with a small Sun angle (Doppler noise ~5x worse), the 
GM of Ceres is not well determined with an uncertainty of about ±0.2 km 3 /s 2 . For the 
same reason, the spherical harmonics solution for the HAMO phase is degraded as shown 
in Fig. 1 1 , where the coefficients are determined through degree 4 and J 2 (unnormalized) 
is determined to about 3% (10 x formal uncertainty). These results are for Doppler tracking 
only with 5.5 h of X-band data per day at 0.5 mm/s accuracy (degraded for the Sun angle). 
The GM and pole determination of the orbit depends on the nonspherical nature of the 
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gravity field. For our simulation a nominal 0.03 J 2 value (unnormalized) is used with the 
Earth gravity field scaled to Ceres used for the remaining coefficients. 

The Ceres LAMO mission phase begins June 22, 2015 and lasts for 30 days in a polar 
circular orbit. The orbit radius is 1180 km resulting in an orbit period of 8.9 hours. The Sun 
angle is more favorable for this part of the mission and so a better accuracy of 0.1 mm/s 
for the Doppler data can be used. With this orbit and tracking accuracy, the gravity field 
can be determined to degree and order 10, or a surface resolution of 140 km (<300 km 
requirement). As in HAMO, these results assume DSN tracking only with 5.5 hours per day. 
The optical imaging landmark tracking scenario has yet to be simulated for Ceres. Flowever, 
similar improvement as with Vesta is expected in the long-wavelength gravity uncertainties. 
The uncertainty in the GM of Ceres reduces in the LAMO orbit to ±0.02 km 3 /s 2 (10 x 
formal), and the Jo uncertainty reduces to 0.4%. Both these results should improve (~10x) 
with optical landmark tracking and if additional Doppler tracking is available. 

The Ceres pole and rotation rate uncertainties expected from the HAMO mission phase 
are given in Table 3. Results are only presented for DSN tracking only, and uncertainties are 
expected to improve similar to Vesta (~5 x for pole and up to 1000 x for the spin rate) once 
landmark tracking observables are incorporated into the analysis. Additional DSN tracking 
via the LGA represents another possible improvement for Ceres. 
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